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ON THE EXCURSIONS OF REFLECTED LOCAL TIME 

PROCESSES AND STOCHASTIC FLUID QUEUES 
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Abstract 

This paper extends previous work by the authors. We consider the local time 
process of a strong Markov process, add negative drift, and reflect it a la 
Skorokhod. The resulting process is used to model a fluid queue. We derive an 
expression for the joint law of the duration of an excursion, the maximum value 
of the process on it, and the time distance between successive excursions. We 
work with a properly constructed stationary version of the process. Examples 
are also given in the paper. 
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1. Introduction 

Consider a stationary strong Markov process X = (Xt ,t G M) , defined on 
some filtered probability space {Vt,T,P,{J-t,t G M)), with values in M+, a.s. 
cadlag paths, and adapted to {Ft). In this paper, the local time L of the 
process X at x = is considered an (J^t)-adapted stationary random measure 
that regenerates jointly with X at every (stopping) time that X hits 0. More 
precisely: 
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(Al) L assigns a nonnegative random variable L{B,uj) to each B G B{M) such 
that L{-,uj) is a Radon measure for each u £ Q. 

(A2) For any a.s. finite (J^t) stopping time T at which Xt = 0, the process 
{{XT+t, L{T,T + t)),t > 0) is independent of J-t- 

We take the broader perspective with regard to the process L and we allow 
for the case that it is a local time of an irregular point (in which case L has 
discontinuous paths) as well as the case that is a sticky point (in which case 
L is absolutely continuous with respect to the Lebesgue measure with density 
c^{Xt = 0) for some c > 0). We refer to [21 Chap. IV] (in particular Corollary 
6), |17t Chap. 6], and [SI §V.3] for further discussion. For each s £ M define the 
inverse local time process with respect to t by 

L;;i := inf{t > : L[s,s + t] > n}, u > 0. (1) 

What is important is that, owing to this definition, the inverse of the cumulative 
local time is a Levy process in the following sense: 

Lemma 1. If L is continuous then for every a.s. finite [Tt)- stopping time T 
such that Xt = 0, the process {L^\,u > 0) is a subordinator with L^^q = 0. 

If L is not continuous, that is to say if is an irregular point for X, then this 
Lemma is taken as an additional requirement to the definition of L. This is 
easily arranged by choosing L to be a modification of the counting process on 
Z, the discrete set of times that X visits 0, so that the inverse is a subordinator. 
To do this, we assign, to each element of Z, an i.i.d., unit-mean exponentially 
distributed weight. Then let the local time on an interval I to be the sum of 
all the weights of the points of Z in /. 

We summarise this as an assumption, in addition to (Al)-(A2) above: 

(A3) If L is discontinuous then we require that for every a.s. finite (J-t) stopping 
time T such that Xj- = 0, the process {L^^-: W>Q) is a subordinator. 

We will also need the following assumption: 
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(A4) The stationary random measure L has finite rate not exceeding 1, i.e. 



EL{0,t) = fit, 

where < < 1. 

Then, as in [19], [16], [22], and [15], we define a stationary process Q = 
iQt,teR) by 

Qt= sup {L{s,t]- {t- s)} , t£R. (2) 

—oo<s<t 

Furthermore, Q is ergodic (its invariant u-field is trivial.) Notice that Q also 
satisfies, pathwise, 

Qt = Qs + L{s,t]-{t-s)- inf (Qs + L{s,r]-{r-s)) AO, 

s<r<t 

= sup (L(r, t]-{t- r)) V (Q, + t] - {t - s)) (3) 

s<r<i 

for all — oo < s < t < oo. It is worth recalling |15j that if we consider ([3|) as a 
fixed point equation for Q then process defined by ([2|) is the unique stationary 
and ergodic solution of ([3|). A typical sample path of Q is depicted in Figure 
[2] below. It consists of isolated excursions away from zero (also called "busy 
periods"), followed by intervals of time at which Q stays at zero (called "idle 
periods"). In this respect, the process Q is thought of as the workload in a 
stochastic fiuid queue. Amongst other things in |16] . [22], and [15] . expressions 
are derived for the marginal distribution of Q and the Laplace transform of the 
duration of a typical idle and busy periods. 

In this paper, we shall derive an expression for the joint law of three random 
variables: the duration of a busy period, the duration of an idle period, and 
the maximum of Q over a busy period. 

It is assumed, throughout, that {Cl,T,P) is endowed with a P-preserving 
measurable fiow : — t- f], t € M, with a measurable inverse 6^^ = 0_j.All 
stationary random processes and measures can be constructed on Q in such a 
way that the fiow commutes with the natural shift, e.g., Qt{9s°oj) = Qj+s(w), 



4 



Takis Konstantopoulos, Andreas E. Kyprianou, Paavo Salminen 



and L{B,9sUj) = L{B + s,uj), for all s,t G M, Borel sets B C R, and uj £ 
The flow will be explicitly used in Section [5] to obtain distributions conditional 
on observing a positive (or a zero) value of Qq. 

2. A closer look at the reflected process 

Consider now any a.s. finite (J-t)-stopping time T, such that Xt = 0. Then 
{Lj^\-i t > 0) is a subordinator starting from zero (owing to Lemma [1] or 
assumption (A3)) with law that does not depend on T. It turns out that 
the process of interest is A^^ = {h.T,t : i > 0}, where 

Ar;t = t - L^l, t > 0. (4) 

Note that, irrespective of T, the process At has the law of the same bounded 
variation, spectrally negative Levy process which is issued from the origin at 
time zero. By (A4), L has rate < 1; hence EAt-^i = 1 — ^ < 0. Since L^^^ 
is a subordinator, it has a well-defined, possibly nonzero, drift. If this drift is 
larger than or equal to unity then — A^^ is a subordinator and, as it will turn 
out, this is a trivial case. 

We therefore assume in the sequel that the drift of L^^ is less than unity or, 
equivalently, that 

(A5) The drift 5/<^ of the process A defined by (jH) is strictly positive. 

Under this assumption, the point is irregular for (— oo,0) for Aj- (this 
follows as a standard results for bounded variation spectrally negative Levy 
processes, see Bertoin [21 Chap. VII].) 

In addition, under (A5), it is clear that the time taken for At to first enter 
(—00,0) is almost surely strictly positive. It will be shown below (Lemma 
[3]) that this implies that the excursions of the process Q, i.e. the busy periods, 
have strictly positive Lebesgue length with probability one. It can be intuitively 
seen, via a geometric argument involving the reflection of the space-time path 
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Figure 1: The construction of the process {Atu, i > 0) and related processes, assuming that 
T = 0. Note that A may have countably many jumps on finite intervals. 

of At about the diagonal (see Figure [1]) , that the time taken for At to first 
enter (— oo, 0) is almost surely equal to the length of the excursion of Q started 
at time T. 

In this light, note also that At cannot creep downwards because it is spec- 
trally negative with paths of bounded variation (cf. Bertoin [21 Chap. VII]). 
Hence the overshoot at first passage of At into (— oo, 0) is almost surely strictly 
positive. It will turn out (Lemma [1]) that this overshoot agrees with the idle 
period following the aforementioned excursion of Q. 

The above analysis implies that, on finite intervals of time, Q has finitely 
many excursions (busy periods) separated by positive-length idle periods. De- 
note by 



■ ■ ■ < g{-l) < g{0) < g{l) < g{2) < ••• 



the beginnings of the idle periods and by 



• • • < d{-l) < d{0) < d{l) < d{2) < ■■■ 
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their ends, see Figure[2l We choose the indexing so that ^(O) < < g{l). Let Ng 
(respectively, Nd) be the point process with points {g{n) : n £ Z} (respectively, 
{d{n) : n G Z}). As Q is a stationary process, Ng and Nd are jointly stationary 
with finite, nonzero, intensity [15] denoted by A (an expression for which is 
given by ()26p and is derived in ^4.31 below). Corresponding to point processes 




d{-l) g(0) d{0) g{l) d{l) 



B I 

Figure 2: The definition of g{n) and d{n). By convention, the origin of time is between g{Q) 
and under the original measure P. Under Pa, the origin of time is at d(0). Under Pg, 

the origin of time is at g(0). The random variable Q* is the maximum deviation from of Q 
within the typical busy period. 

Ng, Nd we have the Palm probabilities Pg, Pd, respectively. Let us consider Q 
under the measure Pd- Then Pd{d{0) = 0) = 1, i.e. the origin of time is placed 
at the beginning of a busy period. By the strong Markov property, the "cycles" 

Q := {Q^ : d{n) < t < d{n + I)} , n € Z, 
are i.i.d. under measure Pd- In particular, the pairs of random variables 

{g{n + 1) - d{n), d{n + 1) - g{n + 1)) , n G Z, 
are i.i.d. under Pd- Consider the triple 

{B,I,Q*)--=(g{l)-d{Q), d{l)-g{l), sup QA , (5) 

\ d(0)<t<g(l) / 

which is a function of Cq. We are primarily interested in the P^^-law of {B, I, Q*) 
Since, under Pd, the origin of time is placed at (i(0), we interpret B,I,Q* as 
the typical busy period, the typical idle period, and the maximum value of Q 
over a typical busy period, respectively. 
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The next lemma is proved in |15] : 

Lemma 2. Let D = mf{t > : Xt = 0} and d = mi{t > : Qt>0}. Then 
d = D a.s. on {Q = 0}. 

We now obtain an alternative expression for B = g{l) — d{0) and / = d{l) — 
g{l) in terms of the inverse local time. 

Lemma 3. We have that 

B = g{l) - d{0) = mf{u > : > u}, (6) 

^ + ^ = ^(l)-^(0)=^,-(0);.(iM(0)- 

Proof. Since d{0) is the end of an idle period (and the beginning of a busy 
period), we have Qd{o)- = 0- Using then expression ([3]) we obtain 

Qt = L[diO),t] -it- dm, d{0) < t < 5(1), 

which gives 

B = g{l) - d{0) = inf{t > : L[d{0),d{0) + 1] = t}. 

Consider now -^^^(o).^, defined by ([1]). By Lemma[2l d{0) is a point of increase of 
the function t h-^ L[(i(0),(i(0) + t]. Hence 5(1) > d{0). Also, when L[d(0), (i(0) + 
t] — t decreases, it does so continuously. Therefore, 

B = inf{t > : L[d{0), d{0) + t] < t}. 
Notice also that, for all t,x > 0, 

L[d(0),d(0) +t]<x ^ t< 
where ^^^(o);^,. = li™e;o -^-^(o);j:_e- I* follows that, 

B = mi{t >0: t< 

= inf{t>0: L^loy^t>t}, 
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by the right continuity of 1 1— t- L^^^^.^. To prove the expression for B + 1, notice 
that L does not charge the interval [^(l), d{l)) because, by definition, Q is zero 
for all t in this interval. □ 

Henceforth it will be convenient to work with the process A = ( ,t > 0) 
where 

Note also that d{0) is an (J"f)-stopping time at which X takes the value and 
hence in our earlier notation At = A^^q).^. 

From the expression ^ , and as discussed in the introduction of Section [2l 
we see that B is simply the first time at which A enters (— oo,0), 

B = inf{t > : Af < 0}, (8) 

which is necessarily strictly positive thanks to the irregularity of for (— oo,0) 
for A. From ^ and ([7]) we see that 

^ = -^d(0);g{l) - ^ = ^diO);B " = "^S, (9) 

i.e. / is, in absolute value, equal to the value of A at the first time it becomes 
negative. Again we recall from the discussion at the beginning of Section [2] that 
A cannot creep downwards and hence / > almost surely. 

Consider now the random variable Q* = sup^(Q)<t<g(]^) Qt- If we define 

Tx := inf{t > : At > x} = mi{t > : At = x} (10) 

we immediately see that 

{Q* <X} = {B< Tx}. (11) 

3. The triple law 

Recall that Pd is the Palm probability with respect to the point process 
{d{n),n G Z}. The function 

H{a,f3,x) = Ed [e-°^-^^l(Q* < x) 
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characterizes the joint law of the triple {B,I,Q*) under P^. Since Pd{d{0) = 
0) = 1, we have that 

Aj = t - Lq.J, with Ao = 0, Prf-a.s. (12) 



Recalling the expressions ([8]), ([9]) and (|TT]) for B, I and Q* , respectively, we 
write 

H{a,f^,x) = Ed [e--«+^^-l(i? < r,)] . (13) 

Since our primary object is the process A defined in (|12|) . and in view of (HD 
and (jl3p . it makes sense to consider the process on its canonical probability 
space and denote its law by Pq. Then 

H{a, P, x) = Eo \e-''^+l^^^ 1{B < r^)j . (14) 

The latter function may now be expressed in terms of so-called scale functions 
for spectrally negative Levy processes. To define the latter, let 

tl;A{9) :=logEoe^^S 0>O, 

be the Laplace exponent of A under Pq. Then the, so-called, g-scale function 
for (A,Po), denoted by W^''\x), satisfies W^'i\x) = for x < and on [0,oo) it 
is the unique continuous (right continuous at the origin) monotone increasing 
function whose Laplace transform is given by 

/•CO -I 

/ e-'-W^'^\x)dx = —- , for/3>$A(g), (15) 

Jo -ipAiO) - q 

where 

$A(g) = sup{9 > : ^a(^) = q} 
is the right inverse of (See for example the discussion in Chapter 9 of |17]). 



Theorem 1. Let A be the process defined by (jl2p . B its first entry time to 
(—00,0) as in ([8]), andTx the first hitting time of{x} as in ([T0|). For a,f3,x > 
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we have 



H{a, 13, x) 



En 



1 l + ia-Ml3))Ioe-l'yw(»\y)dy 
6a e-^^W('^){x) ' ^ ' 



Proof. Let Qt := a{As,s < t) and define, for all /3 > 0, the exponential 
(^t)-martingale 

Mf := e^^^-^-^^^\ t > 0. 

Let, on the canonical space of A, Pq be a probability measure, absolutely 
continuous with respect to Pq on Qt for each t, with Radon-Nikodym derivative 



dl 



Qt 



Notice that A is still a Levy process under Pq with Laplace exponent 

i^lid) = logEy''^ = iJAif] + e)- MP)- 



(17) 



It is straightforward to check from the above formula that, under Pq, A is 
spectrally negative, with bounded variation paths and drift coefficient equal to 



6a. Since on the stopped u-field we have 
substitute 



M^, we may 



in the equation (jl4p for H to obtain 

H{a,f3,x) = Eo [M^e'^^(^)^e-"^1(S < r^)j = L^i'^-'f'^m^ < r^) 
Let 

q:= a- V'a(/3), 

and assume that (7 > 0. It follows from [17[ Thm, 8.1(iii)] that 



H{a,l3,x) = E(^ [e-^^±{B < r,.)] = 4'^(0) - z'^J'\ 



wj.'hx) 



(18) 
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where wj^^^ is the (^-scale function for (A,Pq) and zj^^ is given by 



Z^^\x) = l + q£w^'^\t)dt. 



It is easy to see \17\ Lemma 8.4] that the Laplace transform of W^'^\-) is the 
Laplace transform of W^'^\-) shifted by (3 and this ensures that 

= e-^^VF(°)(x). (19) 

Moreover, since A still has drift coefficient 6a under Pq , \n\ Lemma 8.6] tells us 
that, irrespective of the value of q and /3, W^'^^(O) = 1/5a. Putting the pieces 
together, this gives us the desired expression for a > ip\{/3). However [T71 
Lemma 8.3], since W'^'^^x) is analytic in q, the condition on a can be relaxed 
to a > by using a straightforward analytic extension argument. □ 

In view of dH), dH), ([TT]) . and (fT3|) . we get the following corollary. 

Corollary 1. (Joint law of typical B, I and Q* ■) Assume that (A1)-(A5) hold. 
Then the joint law of the length B of a typical busy period, the length I of a 
typical idle period, and the maximum Q* of Q over the typical busy period is 
expressed by the formula 

E,[e ^ HQ <x)]-l-- e-(^.wi^){x) 

where a, f3,x > 0. 



4. Marginal distributions 

Clearly, formula (|20p can be used to extract more detailed information about 
typical behaviour of Q. Let us first derive the distribution (Laplace transform) 
of the pair {B,I) under the measure Pd- We have 



Edie-'^'^e-^^] =Eo[e-'^^e^^"] = lim F(a,/3, 



x). 
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To derive the limit, let us temporarily assume that q = a — iph{P) > and 
/3 > 0. Consider in the form (fTSj) and use the limiting result 

Z^^\x) q 
lim — = 

from [171 Exercise 8.5] where the function <I>^ is the right inverse of ip^. That 
is to say 

= sup{^ > : = q} 

= snp{e > : yjlie) = a - i^Am 
= sup{9 > : ^a(6I + /3) = a} 

= ^>A(a) - /3. 



This gives 



f211 



To remove the restriction that a > il'lviP) iii pip and replace it instead by just 
a > 0, one may again proceed with an argument involving analytical extension 
taking care to note for the case that a = 'i/'A(/3), 

a-Ml3) y ^^(^a(«)-/3) , 

|a-^A(/3)K0 ^>A(a) -/3 |a-VA(/3)K0 $A(a) - /3 
4.1. Busy period 

Letting /3 = in (|21|) . we find the P^^-law of B. That is to say, 



Ed[e-n = 1 



5a ^'A(a) 

This formula is consistent with the result of [151 Prop. 8] and, moreover, we see 
that the mean duration of the busy period is given by 
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4.2. Idle period 

To find the P^-law of / we need to set a = 0. Recall however from the 
beginning of Section [5] that E^{Ai) < 0. This implies that 'I>a(0) > and 
hence we have 

It follows that the mean idle period is thus equal to 

E4Il = ^^. (24) 

where Va(0+) = ^rf(Ai) < 0. 

4.3. Rates 

A cycle of the process Q is defined as the interval from the beginning of a 
busy period until the beginning of the next busy period. We therefore have 

mean cycle length = Ed[B + I] = ^ 7 ^^(0+) ^ (25) 

Oa^a{0) 

We can express the common rate, A, of Ng and as the inverse of the mean 
cycle length: 

A:=EiV.(0,l) = £A.,(0.1) = ^^ = ^M^^ (26) 

4.4. The maximum over a busy period 

We now derive the P^-distribution of Q* . Letting a = /? = in (fT6]l we 
obtain 

where W{x) = W^'^^x) is defined through its Laplace transform 



POO I 

/ e-^''W(x)dx = — --, for e > 
Jo IpAiO) 



(0). (27) 



An immediate observation is that liuix^oPdiQ* ^ 2;) = 0, since VF(0) = 
limg^Q 9 /'iIj\{9) = 1/(5a- So under P^, the random variable Q* has no atom 
at zero-which is, of course, expected. 
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We now show that Q* has exponential tail under and derive the precise 
asymptotics. To do this, let 

Then (j27p gives that the Laplace transform of x i— )• e~^'^W{x) is i— )■ 1/^a(/3* + 
9). From the final value theorem for Laplace transforms, 

Jnn e-^'^W{.) = hm ^-^^ = 

where we used the fact that ipAiP*) = 0. It follows that 

OA 

as X — >• oo. 



5. Cycle formulae 

We now show how the use of cycle formulae of Palm calculus enable us to 
find (Proposition [1] below) the joint law of the endpoints of an idle period 
conditional on the event that the idle period contains the origin of time. Also 
(Proposition [2] below) we characterise the joint law of the endpoints of a busy 
period, together with the maximum of Q over this busy period, conditional on 
the event that the busy period contains the origin of time. 

It is well-known that if ($7, J-", P) is endowed with a P-preserving flow {6t,t > 
0) (see end of Section [T]) then for any random measure M with finite intensity 
Am, and any point process with finite intensity Aat such that M(B,6tU)) = 
M{B + t,uj), N{B,9tUj) = N{B + t,uj), for t e M, B Borel subset of M, and 
Lo £ and any nonnegative measurable Z : — )• M, we have 

XmEm[Z] = XnEn [ ZoOt M{dt), (28) 

where Pm, Em (respectively, P/v, Ej\[) denotes Palm probability and expecta- 
tion with respect to M (respectively, N), Tq is the first atom of N which is 
< 0, and Tk,Tk+i are any two successive atoms of N. 
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The next result can be found for some special cases in [16] (diffusions), |15j 
(Levy processes), and in |22| the general expression is derived. Here we offer a 
new proof in the general case based on ([28|) . 



Proposition 1. (Joint law of endpoints of idle period.) Assume that (Al)- 
(A5) hold. Then, conditional on Qq = 0, the left end-point, g{0), and right 
end-point, d{0), of the idle period containing t = have joint Laplace transform 
given by 

^ ^ -^^^0+) a-/3V«-$A(0) /3-$A(0)y'' 

for non-negative a and (3 (a ^ P). 

Proof. Let M/ be the restriction of the Lebesgue measure on the idle periods: 
Mi{A)= [ l{Qt = 0)dt, A£B(K). 

J A 

Then E]\jj[Z] = E[Z\Qq = 0] for all nonnegative random variables Z. Apply 
([28]) with M = Mi, N = Na, and Z = e-"'^{o)+/^'i(o): 

fd{0) 

AM.^M.[e-"''(°)+''^(°)] = XEa / ^-adio)oe.+Mo)oe, Mj{dt). 



Jd{-i) 

Here A is the rate of and is given by The rate Xm, is given by 

E,[I] 

- eJbTTy 

Hence 

A 1 5a$a(0) 



Am. Ed[I] -^^(0+)' 
where we used ()24p and ()25p . Now, P^{d{0) = 0) = 1. To compute the integral 
above, note that Mj is zero on the interval {d{ — l),g{0)), and that, for g{0) < 
t < 0, we have d{0)o6^ = —t, and g{Q)°9t = g{0) — t. So the integral above 
equals 

g(a-/3)t-ag(l) « ^ 



9(0) « - /3 
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Combining the above we obtain 

-V'AlO+j a - /3 

Since Ed[e'^9{o)j = Edie-/^^], the result is obtained by using ([23]). □ 

Proposition 2. (Joint law of endpoints of busy period and maximum over it.) 
Assume that (Al)-(A5) hold. Then, conditional on Qq > 0, the left end-point, 
d{0), and right end-point, g{l), of the busy period containing t = 0, together 
with the maximum of Qg for s ranging over this busy period have a joint law 
which is characterised by 



^[e-«9(i)+M0)i(Q* < x) I Qo > 0] 

(29) 



^^{Q) n + aJ^W^^){y)dy _ I + ^ W^P\y)dy \ 



a 

for non-negative a and (3 (a ^ 13) 



Proof. Let Mb be the restriction of the Lebesgue measure on the busy 
periods: 

Mb{A) = [ l{Qt > 0)dt, A G B{R). 

J A 

Then EmAZ] = E[Z\Qo > 0] for all random variables Z > 0. Apply ([28]) : 



\Ed [ e-"^'(^)°'''+'^'^(°)°^'l(Q*o0t < x) Meidt) 

Jd{0) 

XEd / e-"(5(^)-*)-^*l(g* < x) dt 
Jo 



XEd 
A 



a — 13 
A 

a — /3 



HQ* < x) e-'^sW ^ 

a — p 

Ed[e~^3^^^HQ* < ^)] - i^d[e-°^(')l(g* < x)] 
{H{P,0,x)-H{a,0,x)), (30) 
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where H{a,(3,x) is the right-hand side of ([20]). Using ([26]), ([25]) and ([22]), we 
have 



Combining the above we obtain the announced formula. □ 



Proposition[2]yields the next corollary which recovers a result obtained in 
using different methods (for special cases, see |16j and [IS])- Clearly, Corollary 
[2] could also be proved analogously as Proposition [TJ 

Corollary 2. Assume that (Al)-(A5) hold. Then, conditional on Qq > 0, the 
left end-point, d{0), and right end-point, g{l), of the busy period containing 
t = have joint Laplace transform given by 

for non-negative a and [3 (a^ j3). 

Proof. The argument proceeds as in the proof Proposition [2] by omitting 
the factor < x), i.e. by formally replacing x with +oo. The last line of 

(f30|) will give ^j^(i/(/3, 0, oo) — H{a, 0, oo)), where H{a, oo) is given by the 
right-hand side of (HH). □ 

Corollary 3. Assume that (A1)-(A5) hold. Then, conditional on Qq > 0, the 
maximum of Q over the busy period containing t = has distribution 

. * , ^ , .W(x) f^W(y)dy- f^W(x-y)W(y)dy , , 

W\xY 

for X > 0. 

Proof. Letting a, /3 — )■ in (j29l) yields 

OH 

P{Q* <x\Qo>0) = «>a(0) lim -—{a,0,x), 

a^o oa 

where 

H{a,U,x) = 1 , , , , . 
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Next recall that for each x > 0, W^"^ (x) is an entire function in the variable a 
and in particular 

where W*^''~^^\x) is the (A; + l)-th convolution of W (cf. Bertoin [3]). From 
this one easily deduces that 

^T^(")(x)U=o = £ W{y)W{x - y)dy. 
The result now follows from straightforward differentiation. □ 

6. Example: Local time storage from reflected Brownian motion with 

negative drift 

Let X = {Xt,t G M} be a reflected Brownian motion with drift — c < in 
stationary state living on / = [0, oo), and let Pq denote the probability measure 
associated with X when initiated from at time 0. Its local time (at 0) for 
s < t is given by 

L(s,t] := lim— / l[o,,) (X„)(iu. (32) 
Let Q be the stationary process defined as in 

Qt := sup{L(s,t] - {t- s)}. 

s<t 

This particular example of fluid queues was introduced and analysed in [19] 
and further studied in [16] and |15| . 

Recall that EqL{0, 1] = c, and, hence Q is well-defined if and only if < 
c < 1 . Here we make this example more complete by finding the a-scale function 
associated with process At := t — L'j^^,t > 0, where 

L^^ := inf{s : L(0, s] > t}, t> 0. 

is the inverse local time process. As seen from formulae (j3ip and ()16p . the 
a-scale function is the key ingredient needed for computing the distribution of 
the maximum of Q over a busy period and related random variables. 
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To begin with, we recall some basic formulae. When normalising as in ([32 
see [ini pg. 214], [6l pg. 22] it holds that 

Eo(exp{-0L,-i}) = exp|-i /"(I - 6"^")^!= e'^Va j 



where 

1 



GeiO,0) :-- 



VW+c^ - c 



r.2 

is the resolvent kernel (Green kernel) of X at (0,0); see [H pg. 129]. Conse- 
quently, we have 

Eo (expjMJ) = exp |t (^9 - ^^^^ } = exp{t^A(e)}, 

where 



^^(0) — e- + c2 + c, 0>O. 

Recall (cf. p5|) ) that the a-scale function [a > 0) associated with A is 
defined for a; > via 

/ e-S^w^'^\x)dx = —- ; (34) 

Jo V'a(6') - a 

for X < we set W^'^\x) = 0. The 0-scale function is called simply the scale 
function and denoted W. For the next proposition introduce 



2 



oo 



Erfc(x) ■■= ^ I e"* dt, 

and notice that Erfc(O) = 1, Erfc(+oo) = 0, and Erfc(— cxo) = 2. 
Proposition 3. The a-scale function W^"'^ of A is for x >0 given by 



T^H(x) = -(Aie^?^/2Erfc(-Aiv^) 

_A2e^2=^/2Erfc(-A2 7^72)), (35) 

where 

Ai := 1 + V(l-c)2 + 2a, A2 := 1 - ^/{l-cy + 2a (36) 
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In particular, 



W{x) = ^.^ ' . ((2 - c)e(^-"^'^/^Erfc(-(2 - c)^/^) 



{2-cYx/2- 



2(1 - c) 

-ce''"^/^'Eiic{-c^x/2)], (37) 



and W{Q) = 1. 

Proof. From (I34p we have 



OO -| 

g-fe ^r(a)^x^dx = . (38) 

6 -VWT^ + c-a 

To invert this Laplace transform, introduce \ := 29 + . With this notation, 
1 2 



e-VWT^ + c-a A-2\/A + 2(c-a) -c2 

= (39) 

(^/A-Ai)(^/A- A2) 

2/1 1 



Ai — A2 \ ^/X — Ai \/A — A2 , 

where Ai^2 are the roots of the equation — 2z + 2(c — a) — = 0, i.e., as in 

P6p . Next, recaU the fohowing Laplace inversion formula (cf. Erdelyi [71 pg. 
233]) 

£-1 (^J—) = - /3e^^- Erfc(/3V^) (40) 
valid for A — > 0. Since 

/•OO 

Jo 

we obtain using ([1P|) 



00 /"OO 







26^=^2/ ^(")(2y) = ^-1^ (Aie^?^Erfc(-Ai^) - A2e^i ^ Erfc(-A2 V^)) , 



which is formula (j35p . In particular, when a = it holds Ai = 2 — c and A2 = c 
yielding formula ([37]) . □ 

Using the scale function W and the fact 

$a(0) = sup{e > : ^pA{9) = 0} = 2(1 - c) 
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formula ()3ip yields the distribution of the maximum Q* over an observed busy 
period (i.e. over a busy period containing the origin of time). 

Proposition 4. Let < c < 1. The distribution of the maximum Q* over an 
observed busy period of a reflected Brownian motion (with drift —c) local time 
storage is given by 



/o" W{y) {W{x) - W{x - y)) dy 
W^{x) 



P{Q* < x I Qo > 0) = 2(1 - c) 
where the scale function W is given by |37[ ). 

We plot the derivative of ()4ip for c = 1/2 in Figure [3] below. 




(41) 



Figure 3: The density of Q* conditional on {Qo > 0} for the example corresponding to 
Brownian motion with drift — c — —1/2. 



We recall some formulas from 1161 . First 

^/ ed(0)-/35(l)|Q^>0 



Ki - C) 



^2^ + (1 - c)2 + ^2/3 + (1 - c)2 
1 



(V2^ + (l-c)2 + 1 + c)(72/3 + (1 - c)2 + 1 + c) 
--:F{e,l3;l-c), (42) 



and 



E(^e'3(^)-P<''^\Qo = 0^=F{e,/3;c). 



(43) 
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Setting f3 = 6 in the right-hand side of ()42p and (j43p . respectively, we get 

)-dm I ^_ ^ n\ _ 4(1 - c) 



E (e-^(f 
and 



^20 + (1 _ c)2(^2^ + (1 _ c)2 + 1 + c)2 ^ 



(44) 



4c 



V20 + c2(V20 + c2 + 2-c)2' ^^^^ 
Taking the inverse Laplace transform of (j44p (cf. Erdelyi pg. 234]) we obtain 
the density of the length of the busy period g{l) — d{0), given Qq > 0, as 



fg_hiv) = 2(1 - c)e 



-{l-cfv/2 



2v/7r - (1 + c)^;e(^+")'''/2Erfc((l + c)v/V2)) 



Note that the density of the length of the idle period d{0) — g{0), given that 
Qo = 0, is obtained from fg_h{v) by substituting c for 1 — c. In Figure H] we 
have ploted fg^b{v) for three different values of c. We notice also that the mean 




■ c=0.1 c=0.5 c=0.9 



Figure 4: The density of length of the busy period, given that Qo > for three different values 
ofc. 



busy period length has a simple expression: 



E[g{l) - d{0) I Qo > 0] 



2 - c 
(1^ 
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The joint density of d{0) and g{l) is given by 

P(-(i(0) G dx,g{l) edy\Qo>0) = 2(1 - c)e-(i-^)'(^+^)/2 

X (^^2/{7rix + y)) - (1 + c)e(i+^)'("+J')/2Erfc((l + c)^{x + y)/2) 

and, again, the density for {g{0),d{0)) is obtained by substituting c for 1 — c. 

Next we find the density of g{l) (recah that —d{0) is identical in law with 
g{l)) by inverting the Laplace transform (obtained from (j42p by choosing 6 = 
0): 

E (e-'^f I Qo > 



4(1 - c) 



^2/3 + (1 - c)2 + 1 - cj (^V2/5+(l-c)2 + l + 
Letting A := 2/3 + (1 — c)^ we rewrite (f46]) as 

V ' / c VVA + 1-c ^/A + 1 + c/ 



(46) 



From (Ho 



-1 



VA + 1-c \/A + 1 + c, 
= (1 + c) e(i+")'^ Erfc ((1 + c) V^) - (1 - c) e^^-^^)'^ Erfc ((1 - c)VE) . (47) 

Consequently, 

2e"(i-'^)'"/g(i)(2x) = ^ii^(^(l + c)e(^+'^)'"Erfc((l + c)V^) 

- (1 - c) e(i-^)'^- Erfc((l - c) V^)) , (48) 



where /^(i) denotes the density of ^(1) conditioned on {Qq > 0}. From ([18|) we 
obtain 

= ^^"'^'1^^'"'^ ((1 + c) e(^+<^)^^/2 Erfc ((1 + c)^) 

- (1 - c) e^i-'^)'^/^ Erfc ((1 - c) ) . (49) 
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Moreover, the density /^(o) of d{0) conditional on {Qq = 0} is obtained from 
I9|) by substituting c instead of 1 — c : 



fm(^) = ^Y^T" " Erfc ((2 - c)^^ 

- ce'^'^/2 E^fp (^cV^) ) . (50) 

It is striking how similar formulae (|37p and (j50p are. 



Remark 1. The scale function formulae (|35p and (j37p are clearly valid for all 
c > 0. In case c = the process {L{0, t] ; t > 0} is a version of the Brownian 
local time, and the a-scale function VFq"^ of the corresponding process A is 
given by 

Wi''\x) = f(l + VrT2^)e"v^™Erfc f-(l + ^/^) 

2vl + 2a V V / 

- (1 - Vl + 2a) e-^^'^+2° Erfc (-(1 - Vl + 2a) ^/^^ ) . 

In particular, 

Wq{x) = e^^ Erfc(-\/2^). (51) 



In case c = 1 it holds Ai^2 = 1 ± v2a and for a 7^ formula ()37p can be used 
directly. For the 0-scale function we need to take the limit as c — )■ 1 in ([57]) : 

T^i(x) = (1 + x) Erfc(-y^) + -\/ — e-^/2^ (52) 

V vr 

Remark 2. Here we display some formulae for Laplace transforms apparent 
from above and point out a misprint in Erdelyi et al. [7]. 

First, from (j38|) . (|39|) . and (j52p we have the following Laplace inversion 
formula valid for A > 1: 



and this can be "extended" (for a > 0) to 



(53) 



C^^ ( ^ ) = (1 + 2a2x)e"'^'Erfc(-a\/^) + 

V(vA - a)2/ 



(54) 



vr 
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Furthermore, it can be checked that (j54p is vahd for all a < by evaluating 
the Laplace transform of the right-hand side. This can be done term by term 
by using, e.g., Erdelyi et al. O pp. 137, 177] (well-known formulae): 



2 

C fe"'^ Erfc(aV^)) = X'^/^ (A^/^ + a)-\ 



and 



C ( rre"'^ Erfc(aV^) ) = (e"'^' Erfc(aV^) 



2a 

We remark that formula (10) in [7j p. 234: 



(^j-^^^—^^ = 1 - l^fh^ + (1 - 2hx) e*^ (^Erf (\/te) - 1 



(55) 



is not correct since it does not coincide with formula ()54p (for a < 0). Indeed, 
because 

Erf(x) := ^ / e-*'dt, 
V^r Jo 

the right-hand side of (|55p is zero at zero but the right-hand side of ()54p is 1 
at zero. 

7. Further examples 

In the previous example we derived a local time process from a given Markov 
process. However, it is also possible to consider examples where just the local 
time process L, or equivalently the subordinator is specified. Indeed the 
subordinator that will play the role of L^^ in this example has no drift and has 
Levy measure given by 



T{v) T{v) 



26 Takis Konstantopoulos, Andreas E. Kyprianou, Paavo Salminen 

where the constants 99, 7 > and v G (0, 1). Note in particular then that is 
the sum of two independent subordinators, one of which is a compound Poisson 
process with gamma distributed jumps, the other has infinite activity and is of 
the so called tempered-stable type. Clearly 11 also describes the Levy measure 
of —A too. 

According to [9], the process A belongs to the Gaussian Tempered Stable 
Convolution class and moreover, 

7 



\7 + 

for ^ > 0. In particular = 1 and <&a(0) = 93. It is a straightforward exericse 
to show that 

E(Ai) = ^^(0+) = 

7 

and this implies that 

1 ^ 

^ 1 + (/?z>/7 

as required. 

Prom [9] we also know that 

where 

^-d^) -=E r(an + /3) 

n>0 ^ 

is the two parameter Mittag-Leffler function. 

We may now deduce from the theory presented earlier that, for example, 

'^^ l + Y e-(''+^)yy''~^E^^^{Yy'')dy 

and 

7 

1 + f 



PdiQ* > rr) ~ 1 
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